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The recently obtained solutions of Dirac equation in the confining SU(3)-Yang-Mills field 
in Minkowski spacetime are applied to describe the energy spectra of quarkonia (char- 
monium and bottomonium) . The nonrelativistic limit is considered for the relativistic 
effects to be estimated in a self-consistent way and it is shown that the given effects are 
extremely important for both the energy spectra and the confinement mechanism. 

1. Introductory remarks 

Theory of quarkonium ranks high within the hadron physics as the one of central 
sources of information about the quark interaction. Referring for more details to the 
recent up-to-date review, 1 it should be noted here that at present some generally 
accepted relativistic model of quarkonium is absent. The description of quarkonium 
is actually implemented by nonrelativistic manner (on the basis of the Schrodinger 
equation) and then one tries to include relativistic corrections in one or another 
way. Such an inclusion is not single-valued and varies in dependence of the point 
of view for different authors (see, e. g. Ref. 2 and references therein). It would be 
more consistent, to our mind, building a primordially relativistic model so that one 
can then pass on to the nonrelativistic one by the standard limiting transition and, 
thus, to estimate the relativistic effects in a self-consistent way. 

As follows from the main principles of quantum chromodynamics (QCD), the 
suitable relativistic models for description of relativistic bound states of quarkonium 
should consist in considering the solutions of Dirac equation in a SU(3)- Yang-Mills 
field representing gluonic field. The latter should be the so-called confining solution 
of the corresponding Yang-Mills equations and should model the quark confine- 
ment. Such solutions are usually supposed to contain at least one component of 
the mentioned SU(3)-field linear in r, the distance between quarks. Recently in 
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Ref. 3 a number of such solutions has been obtained and the corresponding spec- 
trum of the Dirac equation describing the relativistic bound states in this confining 
SU(3)-Yang-Mills field has been analysed. In this note we should like to apply the 
results of Ref. 3 to description of the charmonium and bottomonium spectra. We 
here solve the inverse problem, i. e. we define the confining gluonic field components 
in the covariant description (SU(3)-connection) for charmonium and bottomonium 
employing the experimental data on the mentioned spectra. 4 As a result, we shall 
not use any nonrclativistic potentials modelling confinement, for example, of the 
harmonic oscillator or funnel type, in particular, because the latter ones do not 
satisfy the Yang-Mills equation while the SU(3)-gluonic field used by us does. In 
our case the approach is relativistic from the very outset and our considerations are 
essentially nonperturbative since we shall not use any expansions in the coupling 
constant g or in any other parameters. 

Further we shall deal with the metric of the flat Minkowski spacetime M that we 
write down (using the ordinary set of local spherical coordinates r,"&,<p for spatial 
part) in the form 



Besides we have \S\ = |det(g^)| = (r 2 sin?9) 2 and < r < oo, 0<i9<7r, 
< ip < 2ir. 

Throughout the paper we employ the system of units with Ti = c — 1, unless 
explicitly stated otherwise. Finally, we shall denote L 2 (F) the set of the modulo 
square integrable complex functions on any manifold F furnished with an integration 
measure while LV,{F) will be the n-fold direct product of L 2 (F) endowed with the 
obvious scalar product. 

2. Preliminary considerations 

2.1. Dirac equation 

To formulate the results of Ref. 3 needed to us here, let us notice that the rela- 
tivistic wave function of quarkonium can be chosen in the form 



with the four-dimensional spinors tpj representing j-th colour component of quarko- 
nium. The corresponding Dirac equation for tp may look as follows 



ds 2 = g^dx» <g> dx v = dt 2 - dr 2 - r 2 {dd 2 + sin 2 tfdtp 2 ) . 



(1) 




Dip = [i ip 



(2) 



where is a mass parameter while the coordinate r makes sense of the distance 
between quarks. 
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From general considerations the explicit form of the operator V in local coordi- 
nates x M on Minkowski manifold can be written as follows 

V = i(Y ® h)E^(d^ ® h - \u»abl a l b ®h- igAJ, a<b, (3) 

where A = A^dx^ , = A^T C is a SU(3)-connection in the (trivial) bundle £ over 
Minkowski spacetime, ^3 is the unit matrix 3x3, the matrices T c form a basis of 
the Lie algebra of SU(3) in 3-dimcnsional space (we consider T a hermitean which 
is acceptable in physics), c = 1, ...,8, ® here means tensorial product of matrices, 
g is a gauge coupling constant. Further, the forms uo a b = u^abdx^ obey the Cartan 
structure equations de a — ui a b A e b with exterior derivative d, while the orthonormal 
basis e a — e^dx^ in cotangent bundle and dual basis E a = E^d^ in tangent bundle 
are connected by the relations e a (Eb) = 5%. At last, matrices 7° represent the 
Clifford algebra of the corresponding quadratic form Qi ;3 

— x-^ X2 x^ in 

C 2 . For this we take the following choice for j a 

f=il ")^=(l ^,6=1,2,3, (4) 



where Cf, denote the ordinary Pauli matrices. It should be noted that, in lorentzian 
case, Greek indices jj,, v, ... are raised and lowered with g^ v of (1) or its inverse g^ u 
and Latin indices a,b,... are raised and lowered by r\ a h = rj ab = diag(l, -1,-1,-1), so 
that ejefo"" = V ab , E£E»g^ = r] ab and so on. 

We can concretize the Dirac equation (2) for tp in the case of metric (1). Namely, 
we can put e° = dt, e 1 = dr, e 2 = rdd, e 3 = rsinfldip and, accordingly, Eq = dt, 
Ei = d r , E 2 = d#/r, E 3 = d v /(rsin#). This entails 

W12 = — d$, W13 = — sin tfdip, W23 = —cosfldip. (5) 

As for the connection A^ in bundle £ then the suitable one should be the confining 
solution of the Yang-Mills equations 

dF = F A A - A A F , (6) 

d*F = *F/\A-Af\*F (7) 

with the exterior differential d = dtdt + d r dr + d$di9 + d v dtp in coordinates t, r, ip 
while the curvature matrix (field strcntgh) for ^-bundle is F = dA + A A A and * 
means the Hodge star operator conforming to metric (1). It is clear that (6) is iden- 
tically satisfied — this is just the Bianchi identity holding true for any connection 
so that it is necessary to solve only the equations (7). 

2.2. SU (3) -confining connection 

In Ref. 3 the black hole physics techniques from Refs. 5 was used to find a set of 
the confining solutions of Eq. (7). For the aims of the given paper we need one 
such a solution of Ref. 3 Let us adduce it here putting T c = A c , where A c are the 
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Gell-Mann matrices (whose explicit form can be found in Refs. 5 ) Then the solution 
in question is the following one 

A* + l=At = - a J- +Al ,-Al + ±=A* = a -l±^- {M+ A 2 ),-^Al = - a -l + A 2 , 

A% + ±=A% = bir+B u -A% + ±=A% = -(b 1 +b 2 )r-(B 1 +B 2 ), —j=A% = b 2 r+B 2 

N> N> . (8) 

with all other A c ^ = 0, where real constants aj, Aj, bj, Bj parametrize the solution, 

and we wrote down the solution in the combinations that are just needed to insert 

into (2). As is not complicated to see, the solution is a configuration describing the 

electric Coulomb-like colour field (components A t ) and the magnetic colour field 

linear in r (components A v ). Also it is easy to check that the given solution satisfy 

the Lorentz gauge condition that can be written in the form div(^4) = 0, where the 

divergence of the Lie algebra valued 1-form A = A^Tcdx^ is defined by the relation 

div(A) = l 9li (VW\!rA») ■ (9) 

v\ s \ 



2.3. Dirac equation spectrum and wave functions 

As was shown in Ref. 3 , after inserting the above confining solution into Eq. (2), 
it admits the solutions of the form 



fa = e^'V" 1 ; J1 /V *7« \ ,3 = 1,2,3 (10) 

with the 2D cigcnspinor 3>j = ^ ^ of the cuclidean Dirac operator on the unit 

sphere S 2 . The explicit form of &j is not needed here and can be found in Refs. 6 
For the purpose of the present paper it is sufficient to know that spinors <&., can be 
subject to the normalization condition 

|$ji| 2 + |$j2| 2 )sintfcfoH? = l, (11) 



i. e., they form an orthonormal basis in L 2 ,^ 2 ). 

The energy spectrum e of quarkonium is given by the relation e = u>i + lo 2 + u>3 
with 

-Aig 2 aibi + (m + on)y/(r>% + 2n x a x + A 2 + g 2 b\{n\ + 2n7aT) 

Ul = Ul (n u li,Xi) = n 2 + 2n iai + A 2 

(12) 

u 2 = u 2 (n 2l l 2 ,\ 2 ) = 
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n\ + 2n 2 a 2 + Ai 



(13) 



, , -A 3 g 2 a 2 b 2 + (n 3 + a 3 ) ^(n 2 + 2n 3 a 3 + A 2 )^ 2 + g 2 b 2 (n 2 + 2n 3 a 3 ) 
u; 3 =u J3 (n 3 ,l 3 ,X 3 ) = n 2 + 2n 3 a 3 + A 2 

(14) 

where Ai = Ai - gB x , A 2 = A 2 + .g(Bi + B 2 ) , A 3 = A 3 - 3B2 , rij = 0, 1, 2, while 
Aj = ±(lj + 1) are the eigenvalues of euclidean Dirac operator on unit sphere with 
lj = 0, 1, 2, ... Besides 



ai = ^A 2 -g 2 a 2 , a 2 = A 2 - g 2 ( ai + a 2 ) 2 , a 3 = \J A 2 - g 2 a\ . (15) 



Further, the radial part of (10), for instance, for ^-component, is given at 
rt\ = by 



F n = dAr^e^ 13 ^ (l - ^ ,F 12 = iC x Br a ^e-^ r [\ + 
while at ni > by 

F n = dAr ai e-^ r 



(16) 



"12 



iCi J Br Ql e~' 3ir 



i-|) ^(nJ + ^L^Cn) 

1 + I) x - 1(ri) -V lL ^ 1(ri) 



(17) 



with the Laguerre polynomials L^(ri), n = 2/3ir, /3i = y/Jj-Q — (wi — 5^4i) 2 + <? 2 &i, 
A = ffbi+/3i, -B = /io + wi-ff^i, Fi = [aifii - gcixiui - g Ax) + gaxbi]B + g 2 a 1 b 1 A, 
Zi = [(Ai — gB\)A + ga\fio)}B + g 2 aibiA. Finally, C\ is determined from the 
normalization condition 



r°° 1 

J (\F 11 \ 2 + \F 12 \ 2 )dr = -. 



(18) 



Analogous relations will hold true for ^2,3 , respectively, by replacing a\ , A\ , bi , B\ , ot\ 
a2,A 2 ,b 2 ,B 2 ,a 3 for ip 3 and a lt A lt h, B lt ai — > -(ai + a 2 ),-(Ai + A 2 ),-(6i + 
frz), -( ^1 + -B2), a 2 for V> 2 so th at (3 2 = ^ - [lo 2 + g{A x + ^2)] 2 + <? 2 (&i + 6 2 ) 2 , 
/?3 = ^/yUg — (W3 — .gA 2 ) 2 + g 2 6 2 - Consequently, we shall gain that Vj G ^ 2 (R 3 ) at 
any ( e 1 and, as a result, the solutions of (10) may describe relativistic bound 
states of quarkonium with the energy spectrum (12)-(14). 

2.4. Nonrelativistic limit 

Before apllying the above relations to a description of charmonium spectrum let 
us adduce the nonrelativistic limits (i.e, at c ^ 00) for the energies of (12)-(14). 
The common case is not needed to us in present paper so we shall restrict ourselves 
to the case of nj = 0, 1 and lj = 0. Expanding u)j in x = we get 



wi(0,0,Ai) = -x— h/ioc 

Ai 
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wi(l, 0, AiJ = -x — + /j, c 



l-^(^) * 2 + 0(* 3 ) 



8 VAi 



which yields at c — > oo (putting % = c = 1 again) 



wi(0,0, Ai) = /i 



1 1 ( gal 
2 I Ai 



,wi(l,0, Ai) = ho 



1 i (ga\ 

8 I Ax 



(19) 



(20) 



Analogously we shall have 



w 2 (0,0, A 2 ) = -no 1 - 
w 3 (0,0, A 3 ) = Ho 



1 ^,a (ai + a 2 ) 
A 2 



x _ l ( ga 2 



2 V A 3 



w 2 (l,0, A 2 ) = -Mo 
,w 3 (l,0, A 3 ) = /io 



1 - 



1 / g(ai + a 2 ] 



, 1 ( 9 a 2 
8 I A 3 



(21) 
(22) 



where, of course, Aj = ±1 and X 2 - = 1. 



3. Relativistic spectrum of charmonium 

Now we can adduce numerical results for constants parametrizing the charmo- 
nium spectrum which are shown in Table 1. 



Table 1: Gauge coupling constant, mass parameter Ho and parameters of the con- 
fining SU(3)-connection for charmonium. 



g 


Ho, GeV 




a 2 


h, GeV 


b 2 , GeV 


B x 


B 2 


0.618631 


3.11409 


0.0102202 


-0.126200 


-3.42927 


-3.91720 


1.99496 


2.10796 



As for parameters A\, 2 of solution (8), only the wave functions depend on them 
while the spectrum does not and within the present paper we consider A\ = A 2 = 0. 

With the constants of Table 1 the present-day levels of charmonium spectrum 
were calculated with the help of (12)-(14) while their nonrelativistic values with 
the aid of (20)-(22) according to the following combinations (we use the notations 
of levels from Rcf. 4 ) 

ri c (lS) : £i = wi(0, 0,-1) + w 2 (0, 0, -1) + w 3 (0, 0, -1) , 

JM1S) : £ 2 =wi(0,0,-l)+W2(0,0,l) + w 3 (0,0,-l) , 
Xco{lP) ■ s 3 = wi(0, 0, -1) + c^ 2 (0, 0, -1) + w 3 (0, 0, 1) , 
XdO-P) ■ £4 = wi(l, 0, -1) + c^ 2 (0, 0, 1) + w 3 (l, 0, -1) , 

7] C {1P) : e 5 = wi(0, 0, -1) + c^ 2 (0, 0, 1) + c^ 3 (0, 0, 1) , 
Xc2(lP) : e e = wi(l, 0, 1) + uj 2 (0, 0, -1) + w 3 (0, 0, -1) , 
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Vc (2S) : e 7 = wi(0, 0, 1) + w 2 (0, 0, -1) + w 3 (l, 0, 1) , 

V>(2S) : £ 8 = Wl (0, 0, -1) + lu 2 (0, 0, -1) + w 3 (l, 0, 1) , 

^(3770) : e 9 = wi(l,0, -1) + w 2 (0, 0, -1) + w 3 (0, 0, 1) , 

^(4040) : £10 = wi(l, 0, 1) + cj 2 (0, 0, -1) + w 3 (0, 0, 1) , 

VK4160) : e n = Wl (l, 0, 1) + w 2 (0, 0, 1) + w 3 (0, 0, 1) , 

^(4415) : e 12 = 0, 1) + W2(0, 0, 1) + w 3 (l, 0, 1) . (23) 

Table 2 contains experimental values of these levels (from Ref. 4 ) and our theo- 
retical relativistic and nonrelativistic ones, and also the contribution of relativistic 
effects in %. 



Table 2: Experimental and theoretical charmonium levels 



e J 


Experimental 
Gev 


Relativistic 
GeV 


Nonrelativistic GeV 
GeV 


Relativistic contribution 

% 


£\ 


2.97980 


2.97980 


3.11255 


4.45516 


£2 


3.09688 


3.09688 


3.11255 


0.506165 


£3 


3.41730 


3.41730 


3.11255 


8.91779 


£4 


3.51053 


3.50702 


3.11967 


11.0450 


£5 


3.52614 


3.53438 


3.11255 


11.9350 


£6 


3.55617 


3.61609 


3.11260 


13.9237 


£7 


3.59400 


3.66485 


3.11967 


14.8758 


£8 


3.68600 


3.71725 


3.11967 


16.0759 


£9 


3.76990 


3.77773 


3.11260 


17.6066 


£10 


4.04000 


4.05359 


3.11260 


23.2138 


£11 


4.16000 


4.17067 


3.11260 


25.3694 


£12 


4.41500 


4.47062 


3.11972 


30.2174 



4. Relativistic spectrum of bottomonium 

Numerical values for the corresponding constants parametrizing the bottomo- 
nium spectrum which are shown in Table 3. With the constants of Table 3 the 



Table 3: Gauge coupling constant, mass parameter ^0 and parameters of the con- 
fining SU(3)-connection for bottomonium. 



9 


Mo, GeV 


ai 


«2 


61, GeV 


b 2 , GeV 


Bi 


B 2 


0.172961 


9.71992 


-0.313881 


-0.188848 


-21.1329 


34.2764 


1.48168 


2.01632 



present-day levels of bottomonium spectrum were calculated with the help of (12)- 
(14) while their nonrelativistic values with the aid of (20)-(22) according to the 
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following combinations (we use the notations of levels from Ref. 4 ) 

T(15) : £i = wi(0, 0, -1) + W2(0, 0, -1) + w 3 (0, 0, -1) , 

Xbo(lP) ■ £2 = wi(0, 0, -1) + lu 2 (0, 0, 1) + w 3 (0, 0, -1) , 

XbiO-P) ■ £3 = wi(0, 0, -1) + co 2 (0, 0, -1) + w 3 (0, 0, 1) , 

Xb2{lP) ■ £4 - Wl(l, 0, 1) + W 2 (l, 0, 1) + w 3 (0, 0, -1) , 

T(2S) : £ 5 - wi(0, 0, -1) + w 2 (l, 0, 1) + co 3 (0, 0, 1) , 

X60 (2P) : £ 6 - wi(l, 0, 1) + w 2 (l, 0, -1) + w 3 (0, 0, 1) , 

Xbi(2P) : £ 7 = wi(l, 0, -1) + w 2 (0, 0, -1) + w 3 (0, 0, 1) , 

Xb2(2P) : £ 8 = wi(0, 0, 1) + w 2 (0, 0, -1) + ws(l, 0, -1) , 

T(35) :£ 9 =wi(l,0,l)+W2(l,0,l)+w 3 (0,0,l) ) 

T(4S) : £io = wi(0, 0, 1) + w 2 (l, 0, -1) + w 3 (l, 0, 1) , 

T(10860):en =wi(0, 0,-1)+ W2(l, 0,1)+ w 3 (l, 0,-1), 

T(11020) : £i2 = wi(0, 0, 1) + w 2 (0, 0, 1) + w 3 (l, 0, 1) . (24) 

Table 4 contains experimental values of these levels (from Ref. 4 ) and our theo- 
retical relativistic and nonrelativistic ones, and also the contribution of relativistic 
effects in %. 



Table 4: Experimental and theoretical bottomonium levels 



e 3 


Experimental 
Gev 


Relativistic 
GeV 


Nonrelativistic 
GeV 


Relativistic contribution 

% 


£l 


9.46037 


9.46037 


9.73716 


2.92574 


£2 


9.8598 


9.8598 


9.73716 


1.24391 


£3 


9.8919 


9.8919 


9.73716 


1.56436 


£4 


9.9132 


9.95542 


9.72034 


2.36137 


£5 


10.02330 


10.0427 


9.70960 


3.31702 


£6 


10.2321 


10.2359 


9.72034 


5.03641 


£7 


10.2552 


10.2521 


9.74790 


4.91837 


£8 


10.2685 


10.2611 


9.74104 


5.06841 


£9 


10.3553 


10.3870 


9.72034 


6.41781 


£10 


10.5800 


10.6315 


9.71349 


8.64370 


£11 


10.8650 


10.8536 


9.71349 


10.5041 


£12 


11.019 


11.0312 


9.74104 


11.6955 
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5. Physical interpretation 

5.1. Role of the colour magnetic field 

The results obtained allow us to draw a number of conclusions. As is seen from 
Tables 2,4, relativistic values are in good agreement with experimental ones while 
nonrelativistic ones are not. The contribution of relativistic effects can amount to 
tens per cent and they cannot be considered as small. The physical reason of it 
is quite clear. Really, we have seen in nonrelativistic limit [see the relations (19)— 
(22)] that parameters &i, 2,-81,2 [see Eq. (8)] of linear interaction between quarks 
vanish under this limit and nonrelativistic spectrum is independent of them and is 
practically getting the pure Coulomb one. As a consequence, the picture of linear 
confinement for quarks should be considered as essentially relativistic one while the 
nonrelativistic limit is only a rather crude approximation. In fact, as follows from 
exact solutions of SU(3)-Yang-Mills equations of (8), the linear interaction between 
quarks is connected with colour magnetic field that dies out in the nonrelativistic 
limit, i.e. for static quarks. Only for the moving rapidly enough quarks the above 
field will appear and generate linear confinement between them. So the spectrum 
will depend on both the static Coulomb colour electric field and the dynamical 
colour magnetic field responsible for the linear confinement for quarks which is just 
confirmed by the relations (12)-(14). For bottomonium the mentioned effects are 
smaller than for charmonium which well corresponds to physics in question - in 
the former case quarks are more massive and relativistic effects should be smaller. 
Also one can say in other words that colour magnetic field splits the primordially 
nonrelativistic spectrum into some fine structure as, for example, the Seemann effect 
does in atomic physics. But, unlike the latter case, for quarkonia these corrections 
are not small. 

5.2. Specification of the wave function form 

One can notice that the form of wave functions (16)-(17) permits to consider, 
for instance, the quantity l//?i to be a characteristic size of quarkonium. Under 
the circumstances, if calculating l//?i in both relativistic and nonrelativistic cases 
then one can obtain, for example, for the charmonium 7^(25) -level that relativistic 
size will be of order r ~ l//?i ~ 0.938854 • 10~ 14 cm while nonrelativistic one 
(i.e. if calculating j3\ at b\ = 0) would be of order r ~ 0.101539 • 10 -11 cm, 
i. e., r /r ~ 108 >> 1. Analogously, e. g., for the bottomonium T(lS)-level 
the conforming quantities will be: r ~ 0.618867 • 10~ 14 cm, r ~ 0.379148 • 10" 13 
cm, r /r <~ 6 > 1. This additionally points out the importance of relativistic effects 
for confinement. 

5.3. Comparison with the potential approach 

One should say a few words concerning the nonrelativistic potential models 
often used in quarkonium theory. 7 The potentials between quarks here are usually 
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modelled by those of harmonic oscillator or of funnel type (i. e. of the form a/r + /3r 
with some constants a and (3). It is clear, however, that from the QCD point of view 
the interaction between quarks should be described by the whole SU(3)-connection 
= A^T C , genuinely relativistic object, the nonrelativistic potential being only 
some component of A c t surviving in the nonrelativistic limit at c — > oo. As is easy to 
show, however, the SU(3)-connection of form A\ = Br 1 , where B is a constant, may 
be solution of the Yang-Mills equations (7) only at 7 = — 1, i. c. in the Coulomb- 
like case. As a result, the potentials employed in nonrelativistic approaches do not 
obey the Yang-Mills equations. The latter ones are essentially relativistic and, as 
we have seen, the components linear in r of the whole A^ arc different from A t and 
related with colour magnetic field vanishing in the nonrelativistic limit. That is why 
the nonrelativistic potential approach seems to be inconsistent. Our approach uses 
only the exact solutions of the Yang-Mills equations as well as in atomic physics 
the interaction among particles (e. g. the electric Coulomb one) is always the exact 
solution of the Maxwell equations (the particular case of the Yang-Mills equations). 

6. Concluding remarks 

The calculations of the present paper can be extended. Indeed we have the 
explicit form (10) for the relativistic wave functions of quarkonium that may be 
applied to analysis of the quarkonium radiative decays and electromagnetic tran- 
sitions. Unfortunately in connection with death of E. Choban the work in this 
direction is not yet completed but the first author hopes to discuss the mentioned 
questions elsewhere. 
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To the memory of E. A. Choban 
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